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Ferromagnetic Properties of ZrZn 2 
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The low Curie temperature {Tc ~ 28 K) and small ordered moment (Mo w 0.17/iBf.u. -1 ) of 
ZrZn2 make it one of the few examples of a weak itinerant ferromagnet. We report results of 
susceptibility, magnetization, resistivity and specific heat measurements made on high-quality single 
crystals of ZrZn2. From magnetization scaling in the vicinity of Tc (0.001 < |T — Tc\/Tc < 0.08), 
we obtain the critical exponents /3 = 0.52 ±0.05 and 8 = 3.20 ±0.08, and T c = 27.50 ± 0.05 K. Low- 
temperature magnetization measurements show that the easy axis is [111]. Resistivity measurements 
reveal an anomaly at Tc and a non- Fermi liquid temperature dependence p(T) — po + AT n , where 
n — 1.67 ± 0.02, for 1 < T < 14 K. The specific heat measurements show a mean-field-like anomaly 
at Tc- We compare our results to various theoretical models. 

PACS numbers: 75.50.Cc, 74.70.Ad, 74.25.Fy, 74.70.-b 



INTRODUCTION 



Ferromagnetism in the cubic Laves compound ZrZn2 
was discovered by Matthias and Bozorth— in 1958. Its 
occurrence in ZrZn2 is unusual because neither elemen- 
tal Zr nor Zn is magnetically ordered. The low Curie 
temperature and small ordered moment of ZrZn 2 make 
it one of the few examples of a small-moment or weak 
itinerant ferromagnet. ZrZn 2 was initially considered to 
be a candidate for Stoner theory 3 . However a quanti- 
tative comparison of Stoner theory with experiment sug- 
gests that spin fluctuation effects are important^. In par- 
ticular, the Curie temperature is strongly renormalized 
downwards from the Stoner value estimated from band 
structure parameters^. 

ZrZn 2 crystallizes in the C15 cubic Laves structure 
shown in Fig.l, with lattice constant a = 7.393 A. The Zr 
atoms form a tetrahedrally coordinated diamond struc- 
ture and the magnetic properties of the compound derive 
from the Zr 4d orbitals, which have a significant direct 
overlap leading to the magnetic moment being spread 
out over the network shown by the thick lines in Fig. ^ 
ZrZn 2 is strongly unsaturated: an applied field of 5.7 T 
results in a 50% increase in the ordered moment. In 
contrast, strong ferromagnets such as Fe and Ni show a 
negligible increase of the ordered moment with field af- 
ter a single domain is formed. The unsaturated behavior 
of ZrZn 2 indicates a large longitudinal susceptibility and 
the presence of longitudinal spin fluctuations^. Further 
evidence for the existence of strong spin fluctuations in 
ZrZn 2 is provided by the remarkably large effective mass 
of its quasiparticlesi. The low-temperature specific heat 
cocfficien1j£ and the calculated DOS at the Fermi Leveli 
imply an average mass enhancement 1 ± A ss 5 at zero 
applied magnetic field^. This is the largest known aver- 
age mass enhancement for a ci-band metal, and is even 
slightly larger than that of the strongly-correlated oxide 
system Sr 2 Ru0 4 for which 1 + A « 3.6 [Refs. EHJ. 




FIG. 1: (Color online) The C15 structure of ZrZn2. Zr atoms 
occupy the sites of the larger spheres. 



It has been known for some time that the ferromag- 
netism of ZrZn 2 is extremely sensitive to pressure^. Re- 
cent experiments on the samples studied here^ have 
shown that a pressure of p c — 16.5 kbar causes the fer- 
romagnetism to disappear with a first order transition. 
Thus we may also view ZrZn 2 as being close to a quan- 
tum critical point (QCP)iiiii. In view of the strong lon- 
gitudinal fluctuations present in ZrZn 2 and its proximity 
to a QCP, it has been proposed as a candidate for mag- 
netically mediated superconductivityi£*i§iiLi£. However, 
as discussed elsewhere^, we find no evidence for bulk su- 
perconductivity at ambient pressure in these samples. 

There are surprisingly few measurements of the fun- 
damental properties of high-quality samples (RRR>100) 
of ZrZn 2 in the literature. In this paper, we present a 
study of the transport and thermodynamic properties 
of ZrZn 2 . In particular, we have measured the magne- 
tocrystalline anisotropy, the magnetization isotherms for 
2 < T < 40 K and < B < 5T; the resistivity from 
low temperatures through the Curie temperature; and 
the specific heat capacity for 0.3 < T < 40 K. These 
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properties contain information about the quasiparticles 
and magnetic interactions in ZrZn 2 . For example, the 
temperature dependence of the resistivity at low temper- 
atures gives information about the fundamental excita- 
tions that scatter quasiparticles, and the comparison of 
specific heat and magnetization data gives information 
about the importance of spin fluctuations in ZrZn 2 . 



II. EXPERIMENTAL DETAILS 
A. Crystal Growth and Sample Quality 

ZrZn 2 melts congruently at 1180°C [Refs. • At 

this temperature zinc has a vapor pressure of about 
10 bars and is an aggressive flux. Thus we chose to grow 
ZrZn 2 by a directional cooling technique 2 ^. Stoichiomet- 
ric quantities of high-purity zone-refined Zr (99.99%, Ma- 
terials Research MARZ grade) and Zn (99.9999%, Metal 
Crystals) were loaded into a Y2O3 crucible. The total 
charge used was 4.2 g. The crucible was sealed inside a 
tantalum bomb which was closed by electron beam weld- 
ing under vacuum. The assembly was heated to 1210°C 
and then cooled through the melting point at 2°C hr _1 . 
The ingot was then annealed by cooling to 500° C over 
a period of 72 hr. This method produced, on occasions, 
single crystals of volumes up to approximately 0.4 cm 3 . 
Single crystals produced in this way had residual resistiv- 
ities as low as po— 0.53 pflcra corresponding to a residual 
resistance ratio RRR = p(293K)/ ( o(T -> 0)=105. With 
the exception of Ref. l23l previous reports by other groups 
of the fundamental transport and thermodynamic prop- 
erties of ZrZn 2 have been carried out on samples with 
RRR < 45. 

The residual resistivity po and the Dingle temperature 
determined from de Haas- van Alphen measurements may 
be used to estimate the quasiparticle mean free path 
I due to impurity scattering. For a crystal with cubic 
symmetry 2 *, 

1 e2rl f A Q 1 ilAQ m 
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From band structure calculations^, we estimate the sum 
of the Fermi surface areas to be Sf — 1.9 x 10 21 m~ 2 . 
Hence ^trans = 1350 A. A second estimate of the quasipar- 
ticle mean free path can be made from the de Haas-van 
Alphen effect; values are in the range ^dHvA = 1500- 
2800 A depending on Fermi surface orbit, in approxi- 
mate agreement with I t rans . In general, one expects that 
TdHvA < Ttrans since Tt ra ns is weighted towards large mo- 
mentum changes^, whereas TdHvA weights all scattering 
equally. However, the opposite situation t dHvA > i~ trans 
may also arise in an inhomogeneous sample, since the ex- 
ponential scattering rate dependence of the Dingle fac- 
tor Rf> = exp(— imib/eBT) causes ^dHvA to be strongly 
weighted towards high quality regions of the sample. 
Given these considerations, the two mean free path es- 



timates ^dHvA and ^trans are as consistent as can be ex- 
pected. 

B. Transport and Thermodynamic Measurements 

Resistivity measurements were made using a standard 
a.c. technique using a Brookdeal 9433 low-noise trans- 
former and SR850 digital lock-in amplifier with a mea- 
suring frequency f—2 Hz. Sample contacts were made 
with Dupont 4929 conducting Ag/epoxy. Measurements 
of a.c. susceptibility were made by a standard technique 
in which the sample was mounted inside a small coil of 
approximately 2500 turns. The system was calibrated 
using the superconducting transition of an indium sam- 
ple of similar size to the ZrZn 2 sample. Magnetization 
measurements were made using a commercial Quantum 
Design MPMS-XL SQUID magnetometer. 

Heat capacity measurements were made both by a 
long-pulse method and an a.c. method 26 . In the long- 
pulse technique 2 ^ the sample was mounted on a silicon 
platform connected to a temperature-controlled stage by 
a thin copper wire. In the a.c. technique the sample was 
mounted on a flattened 12 pva chromel-constantan ther- 
mocouple and heated optically. The long pulse technique 
allows an accurate determination of the absolute value of 
C(T) but has relatively low resolution. The a.c. tech- 
nique has much higher resolution (~ 0.1%), and is thus 
ideally suited to resolving the small anomaly at Tc ■ How- 
ever, it does not allow an accurate determination of the 
absolute value of C(T) due to the poorly defined addenda 
contribution. We estimated the addenda by measuring 
a Cu sample with similar C(T); the a.c. data with the 
estimated addenda subtracted were multiplied by a scale 
factor to match the long-pulse data at a single tempera- 
ture. 



III. RESULTS 

A. Magnetization Isotherms and Scaling near the 
Ferromagnetic Transition 

1. Low Temperature Isotherms 

Since the discovery of ZrZn 2 there have been 
many studies of the magnetic properties. Previous 
work 1 ! 12 ! 18 . 23 ! 28 ! 29 ! 30 . 31 ! 32 ! 33 ! 34 has shown that the Curie 
temperature and the ordered moment, are strongly de- 
pendent on sample quality and composition 2 ^. Thus the 
purpose of the present magnetization measurements is 
to attempt to characterize the magnetic properties in 
the clean limit. The high-quality of our samples is in- 
dicated by their low residual resistivity (RRR=105) and 
long mean free path, which has allowed much of the Fermi 
surface to be observed by the dHvA effect. 

Although many plots of the magnetization isotherms 
for ZrZn 2 have been reported 2 ^* 2 ^ 9 ^ improved instru- 
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FIG. 2: Arrott plot of magnetization isotherms in the tem- 
perature range 1.8 < T < 40 K. In the Stoner theory of a weak 
itinerant ferromagnet, M 2 is a linear function of H/M. Inset: 
the derivative d(M 2 )/d{n H/M) of the 1.8 K and 10 K Arrott 
curves; the 1.8K curve is offset by +5x 10" 3 ^| f.u. -3 T" 1 . 
Dotted lines are guides to the eye. Note the feature in the 
1.8 K curve, indicated by an arrow. 



mentation and higher sample quality have recently al- 
lowed fine structure to be observed. For a weak ferromag- 
net, the magnetization isotherms are generally expected 
to obey the Stoner form in which M 2 is a linear function 
of H/M. Fig. [3 is an Arrott plot of our magnetization 
isotherms. The near linearity of the isotherms confirms 
that they are indeed of approximately the expected form. 
However, the isotherm at T — 1.8 K shows an intriguing 
feature at ^lqH w 2.4 T, indicated by an arrow in Fig. El 
Another more pronounced feature is directly visible in 
the T = 1.75 K M(H) isotherm plotted in Ref. [3 at 
a field w 6 T, outside the range of our present measure- 
ments. Structure in the electronic DOS close to the Fermi 
level is expected to have a profound influence on M(T, B) 
in itinerant ferromagnets^*^, and calculations? suggest 
that just such structure is present in the DOS of ZrZii2, 
the Fermi level lying between two sharp peaks in the ma- 
jority spin DOS. This is an interesting connection that 
could be explored in future work. 
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FIG. 3: Low field magnetization isotherms near the ferromag- 
netic critical temperature of ZrZn2 (near-critical isotherm at 
T = 27.6 K shown by filled circles). 
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FIG. 4: Scaling plot of magnetization isotherms of ZrZn2 
near Tc- The data shown here are the same as those in Fig.|3] 
except that we have excluded the region fj,oH < 0.005 T to 
avoid systematic errors due to multiple ferromagnetic domains 
inside the sample (see text). The correct choice of Tc, 13 and 
8 causes the scaled data to collapse onto a single curve for 
T < Tc and another for T > Tc- The values that allow the 
best fit to the scaled data of a 7 th order polynomial in M* are 
shown in Table|I] The solid lines show the best fit polynomials 
for T < Tc and T > Tc- The inset shows the same data on 
log- log axes. 
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2. Scaling near Tc 

At temperatures approaching the Curie temperature 
of a ferromagnet the magnetization isotherms M(H) be- 
come highly nonlinear making the identification of Tc dif- 
ficult. An accurate determination of Tc can be achieved 
by a scaling analysis of the magnetization close to the 
transition; we have therefore measured magnetization 
isotherms at a set of temperatures near Tc allowing us to 
determine Tc and the magnetization scaling exponents. 

Close to the critical temperature of a ferromagnetic 
transition, the magnetization M(H,T) can be expressed 
as M*(h*) where M* = Mt~ !i and h* = (H/T)^! 36 are 
appropriately scaled quantities^. Here t= \1 — T/Tc\ is 
the reduced temperature and /?, 5 are the critical expo- 
nents. The experimental determination of (3 and S allows 
remarkably universal conclusions to be drawn about the 
physical model which underlies the phase transition. We 
have measured the magnetization M{H) in applied fields 
< hqH < 0.2 T at closely spaced temperatures in the 
range 0.1 < \T-T C \ < 2.2 K. Fig.EJshows the raw data. 
In order to determine Tc and the scaling exponents, M* 
and h* were calculated from the experimental data set for 
all values of /3, 6 and Tc within a certain volume of the 
3D parameter space encompassing both the Heisenberg 
and mean field models. For each combination of values, 
the polynomial in M* which best fit the scaled exper- 
imental values of h* was found, and its goodness-of-fit 
(x 2 ) value calculated. The correct parameter values are 
taken as those for which the best fit polynomial has the 
lowest x 2 - Only data points for which noH app > 0.005 T 
are included in the scaling analysis in order to ensure 
the formation of a single ferromagnetic domain inside 
the sample. The data have been corrected for the small 
demagnetizing factor of this sample, D = 0.060, but we 
note that the final results are barely affected by this. 
We find T c = (27.50 ± 0.05) K, (3 = 0.52 ± 0.05 and 
8 = 3.20 ± 0.08 (see Table 0. The values of (3 and 5 
are close to those obtained on lower quality samples with 
smaller Tc's2£. 

In order to interpret our scaling results, we need to 
know whether we are in the critical region. The Ginzburg 
criterion^, 

allows us to estimate the extent of the critical re- 
gion. From neutron scattering measurements^ of the 
wavevector-dependent magnetic susceptibility x(l) an d 
low temperature magnetization measurements we esti- 
mate the magnetic correlation length £o = 33 A. The 
specific heat jump is AC « 155 mJ K _1 mol -1 (see sec- 
tion llllDjl . Hence AT G = 0.4 mK. Thus our data are col- 
lected outside the critical region where mean field behav- 
ior is expected. Table [I] shows that our results are indeed 
in agreement with the expected mean field results for a 
3D ferromagnet. Finally, is it worth commenting that de- 
termining the Curie temperature T c from critical scaling, 




FIG. 5: Arrott plot of magnetization isotherms near Tc for 
ZrZn2. The isotherms were measured at 0.4 K intervals and 
at 0.1K intervals close to Tc; the bold line is the critical 
isotherm at T = 27.5 K. In the Stoner theory of a weak itin- 
erant ferromagnet, M 2 is a linear function of H/M 



TABLE I: Theoretical critical exponents for various models 
of phase transition and the experimental values found in this 
work. T c was found to be 27.50±0.05 K. 
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3.20 ± 0.08 



as described here, yields a lower value than that obtained 
by taking the peak in (dM/dB) with T. This may ex- 
plain the slightly higher values reported elsewhere^ for 
similar samples. 

B. Magnetocrystalline Anisotropy 

The exchange interaction in an itinerant electron sys- 
tem is often modelled as isotropic, depending only on the 
relative orientation of electron spins. However, the crys- 
tal field also enters the free energy via spin-orbit coupling, 
causing anisotropy of magnetic properties. One moti- 
vation for studying the magnetocrystalline anisotropy is 
that it has implications for the symmetry of the supercon- 
ducting order parameter in the spontaneous mixed state 
of the putative superconducting phase of ZrZn2 [Ref . l4fj| . 

Fig. [7] shows a low-field hysteresis loop measured for 
applied fields in the range — 0.05T< noH app < 0.05 T. 
The main panel shows that a field of approximately 5 mT 
is sufficient to create a single ferromagnetic domain. The 
inset reveals the low coercive field of the present samples 
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FIG. 6: The temperature dependence of the magnetization 
(filled circles) and the inverse susceptibility (open squares) as 
determined from Arrott plots. 
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FIG. 7: Magnetic hysteresis curve for ZrZn2 measured at 
T — 5K. The inset shows the central region magnified. Note 
the small coercive field H coal ~ 0.05 mT. 



-ffcocr ~ 0.05 mT and justifies our method (see below) for 
obtaining the magnetic anisotropy constants. 

We have determined the ferromagnetic anisotropy of 
ZrZn2 from M(H) isotherms at T = 5K. A disc-shaped 
sample (diameter 2R = 2.72 mm, thickness t — 0.47 mm) 
was spark-cut so that the plane of the disc was (110). 
This geometry allows access to the three major cubic 
symmetry directions [100], [110] and [111] with the mag- 
netic field in the plane of the disc so that the demagnetiz- 
ing field is the same in each orientation. For this sample, 
the volume averaged demagnetizing factor is D rj 0.146 
[Ref . l4l| . which for M = 0.16 /xb f.u. -1 corresponds to an 
average demagnetizing field /io-Ho ~ 0.037 T inside the 
sample. In this section we denote the measured magneti- 
zation by M|| to emphasize that the SQUID measurement 
is only sensitive to the component of M parallel to H app , 
which is vital for extracting the anisotropy constants by 
the thermodynamic method described later. Fig.[S]shows 
Mm (.Hint) measured with H app parallel to each of the 
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FIG. 8: Magnetization curves of a ZrZn2 single crystal mea- 
sured at T = 5K with field applied along [100] (dotted line), 
[110] (solid line) and [111] (dashed line). The ferromagnetic 
easy axis is seen to be [111] at this temperature. Upper inset: 
plot of Eq.(2) using values of K\ and K2 determined from 
Mm (.Hint) as described in the text. Lower inset: M«(Hi n t) up 
to 5 T for the same sample. 



three symmetry directions; Hint is the internal field that, 
in a simple scalar notation, is approximately related to 
the applied field H app by Hi„ t = H app - DM. As H in t 
increases from zero, a single domain is rapidly formed. 
The slower increase in Mii(ifint) for ^o^int > 0.003 T 
with H app || [001] or [110], arises from the rotation of 
M away from the easy axis towards Hi„t as the ratio of 
the interaction term M • Hi nt to the anisotropy energy 
increases. Fig. [H] shows that saturation is reached most 
rapidly with H app || [111], demonstrating that [111] is 
the ferromagnetic easy axis at this temperature. 

The principal cause of the rounding of the M||(Hi nt ) 
curve, even for H | [111], is likely to be the inhomo- 
geneity of the demagnetizing field. The only practical 
sample shapes that give a uniform demagnetizing field 
are ellipsoids of revolution. We have calculated the vol- 
ume distribution of the demagnetizing field for a perfect 
cylinder of the required aspect ratio assuming that the 
magnetization is uniform; the characteristic width of the 
distribution is SD ~ 0.07 and, for example, ~ 12% of 
the sample volume has a local demagnetizing field that 
exceeds the volume average value by a factor > 1.5. An- 
gular misalignments of up to 2° may also contribute. 

The first two terms in the standard expansion of 
the anisotropic free energy of a ferromagnet with cubic 
symmetry^ can be expressed, for the present geometry, 
as 

F{6) = sin 2 9 (cos 2 9+j sin 2 6)K 1 + - (sin 4 9 cos 2 9)K 2 

(3) 
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where 6 is the angle within the (110) plane between M 
and [100] and K i 2 are the first and second anisotropy 
constants. It is possible, in principle, to determine the 
anisotropy constants by assuming that |M| is fixed in the 
'approach to saturation' region of the Mn(Hi nt ) curved. 
In this picture, the direction of M (and hence the mag- 
nitude of My) for any given Hj nt is determined from 
the condition that the torque acting on M due to the 
anisotropy exactly balances the torque due to the inter- 
action of M with H; nt . However we choose an alterna- 
tive thermodynamic method^ to extract anisotropy con- 
stants from M(H) which is less reliant on assumptions 
about the demagnetizing field. We calculate the total 
magnetic work done in order to bring the sample to sat- 
uration with H app applied parallel to a certain crystal di- 
rection. The work required to produce an infinitesimal in- 
crease, dm, in the sample moment is dW = /ioH app .dm. 
Since the irreversibility of M(H) is negligible in ZrZii2, 
we can relate the magnetic work done to a change in the 
appropriate thermodynamic function of state, namely the 
Gibbs free energy, and we obtain AG = J p; H app ■ dm. 
The anisotropy energy F(0) is therefore given (apart from 
an irrelevant constant) by the area under the H apv (M\\) 
curve. From Eq. (2) it follows that K\ = 4(i r j 110 ] —^[001]) 
and K 2 = 9(F [001] + 3F [m] - 4F [110] ). Evaluating 
the areas corresponding to the magnetic work, we find 
^[110] --F[ui] = 1-1 x lO^eVf.u.- 1 and F [m - F [in] = 
1.6 x 10~ 8 eVf.u. -1 . The cubic anisotropy constants are 
therefore K x = -2.0 x 10~ 8 eVf.u.~\ K 2 = -2.6 x 
10- 7 eVf.u. _1 . The magnitude of the ferromagnetic 
anisotropy is therefore a factor ~ 3 smaller than that 
found in the weak itinerant ferromagnet M3AI [Ref. 
(lerg cm~ 3 = 3.15 x 10~ n eVf.u. -1 in ZrZn 2 ), where 
the relative sizes of K\ and K 2 are opposite to those in 
ZrZii2, and almost two orders of magnitude smaller than 
in Ni [Ref. lij. 

The magnetocrystalline anisotropy of ferromagnetic 
metals can in principle be obtained from band-structure 
calculations of the total energy in which the effect of spin- 
orbit coupling is included^*^. It would be interesting to 
see whether such calculations correctly predict the mag- 
netic anisotropy of ZrZii2. 



C. Resistivity 

Many systems close to a quantum critical point have 
been shown to exhibit a so-called non-Fermi liquid tem- 
perature dependence of the resistivity. Non-Fermi liq- 
uid is generally accepted to mean that the low temper- 
ature exponent n in the equation p(T) = po + AT n is 
not equal to 2. Notable examples of systems exhibit- 
ing non-Fermi liquid power laws in the resistivity include 
high-temperature superconductors (n ~ 1 close to opti- 
mal doping), h eavy Fermion antiferromagnets, such as 
CePd2Si2 [Ref. 47] (n w 1.2 at high pressure close to 
quantum criticality) , and the helical ferromagnet MnSi 
[Ref. 0] (n — > I as T — > 0, in the paramagnetic state 



close to the pressure-induced ferromagnet-paramagnet 
QPT). The low moment ferromagnets N13AI and YN13 
[M(0,0) w 0.075 and 0.04 /i B f .u." 1 respectively] exhibit 
a non-Fermi liquid resistivity over a wider range of pres- 
sures about the QCP and, even at ambient pressure, a 

T 165 dependence is observed over a small temperature 

49 

range . 

Fig. shows the temperature dependence of the resis- 
tivity for ZrZn2. On cooling the sample through Tc, we 
observe a slight 'kink' in the resistivity near the Curie 
temperature. Similar behavior is observed in other itin- 
erant ferromagnets such as nickel, although the behavior 
in ZrZn2 is not as pronounced as that in stronger fer- 
romagnets. The resistivity anomaly associated with the 
ferromagnetic transition can be emphasized if a back- 
ground variation of the form po + AT 2 is subtracted from 
the data; the result is shown in the lower inset. It has 
been argued that the dominant magnetic contribution to 
the resistivity is due to short range spin fluctuations and 
therefore dp/dT should vary like the magnetic specific 
heat near the critical poml^. To test this hypothesis we 
have plotted dp/dT in the upper inset of Fig. |5J There 
is good qualitative agreement with the measured specific 
heat anomaly (see inset to Fig. I10|) . 

At lower temperatures, in the ferromagnetic state, the 
resistivity has a well defined non- Fermi liquid behavior. 
Fitting data over the temperature range 1 < T < 14 K, 
we find n = 1.67 ± 0.02. The inset to Fig. shows 
that when p is plotted against T 5//3 , a linear behavior 
is obtained up to about 15 K. The existence of a T 5 / 3 
power law very close to the ferromagnetic QCP has al- 
ready been established in high-pressure measurements on 
ZrZn2 [Ref. |32( but our result shows that the unusual 
power law applies well away from the QCP. Early data 
on a polycrystalline sample 30 showed a T 5 / 3 temperature 
dependence in the range 15 < T < 50 K but notably not 
at lower temperatures. 

It has been pointed out by a number of workers that a 
p oc T 5 / 3 behavior may be understood in terms of scat- 
tering of quasiparticles by spin-fluctuations^i. A T 5 / 3 
dependence is expected near a 3D ferromagnetic quan- 
tum critical point where the collective spin excitations 
are expected to be overdamped. Thus it appears that 
ZrZii2 is sufficiently close to quantum criticality to ob- 
serve the 5/3 non- Fermi liquid behavior. 

It has been reported that ZrZn2 displays superconduc- 
tivity below about 0.6 K. We find no traces of supercon- 
ductivity in our resistivity curves at ambient pressure in 
samples etched in a HF/HNO3 solutionis. 



D. Heat Capacity 

Fig. ^3 shows specific heat results for T < 40 K on 
single crystal samples that were cut from a region of the 
ingot next to that used for resistivity measurements. As 
with previous measurements 8 a large linear contribution 
to the specific heat C = 7T was observed at low tern- 
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FIG. 9: Upper panel: raw p(T) results for 0.5 < T < 32 K; 
upper inset: dp/dT close to Tc', lower inset: p(T) results close 
to Tc after subtraction of a smoothly varying background 
(see text). Lower panel: p(T) plotted against T 5//3 ; inset: the 
same data over a limited T range, demonstrating the non- 
Fermi liquid power law p(T) = p + AT 5/3 . 



peratures, with 7 = 45 mJ K~ 2 mol -1 . At higher tem- 
peratures the high-sensitivity of a.c. specific heat mea- 
surements allows the small ferromagnetic anomaly to be 
observed. 

The anomaly in C(T) at the ferromagnetic transition 
is only 2.5% of the total in ZrZn2; a reliable separation 
of the magnetic heat capacity C m therefore means that 
both the electronic and phonon components must be de- 
termined precisely, and this is not a trivial task. In order 
to emphasize the ferromagnetic anomaly, in the inset to 
Fig. EH we have plotted the heat capacity after subtract- 
ing a smoothly varying estimate of the non-magnetic heat 
capacity of the form Cbkgd(T') = jT + Cn c byc(T '/©); here 
CDebyo is the Debye heat capacity function and is a De- 
bye temperature. We set 7 = 38 mJ K -2 mol -1 , close to 
that observed at low temperatures, in order to obtain a 




FIG. 10: Specific heat C(T) of single crystal ZrZn2: sym- 
bols show data obtained using a pulse-relaxation technique; 
solid line shows the results obtained by a high sensitivity a.c. 
method. Inset: solid line shows a.c. data close to Tc after 
subtracting a smoothly varying background (see text); dashed 
line shows C m (T) calculated by Stoner theory, multiplied by 
a scale factor 0.5. 



fiat C m (T) immediately above Tc; the only other free 
parameter was allowed to vary to give the best fit to 
the data in a small temperature range 29 < T < 33 K 
above Tc, giving a value Q = 284 K. We emphasize that 
this estimate of the non-magnetic heat capacity may not 
be correct in detail (e.g. there may be significant fluctu- 
ation heat capacity in the range of the fit and we have 
not taken into account the true phonon DOS), but our 
quantitative discussion will be confined to the height of 
the discontinuity at Tc which is hardly affected by the 
choice of background. 

The jump in C at Tc measured here is about a factor 
of 2 larger than that measured previously^ on a sam- 
ple with Tc = 10 K. The shape of the anomaly, shown 
in the inset to Fig. ^3 broadly resembles that expected 
for a mean-field second order transition. However, closer 
inspection shows that there is curvature in our estimate 
of C m (T), both above and below the jump at Tc, up to 
1.5 K from Tc- This almost certainly results from ther- 
mal fluctuations, although our uncertainty in the back- 
ground precludes a detailed analysis. There is also evi- 
dence of some rounding of the anomaly from sample in- 
homogeneity on a smaller temperature scale ~ 0.2 K. 

We now compare our results with Landau and spin 
fluctuation theories. In the Landau approach, the free 
energy is written in the form 



F = -a(T)M 2 



1, 



-6M 4 - MB. 



(4) 
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First, we relate the parameters of the theory, namely the 
Landau coefficients a(T) and b, to experimental quan- 
tities determined near Tc- In mean field theory one 
usually assumes that the coefficient a(T) varies linearly 
with temperature near Tc so that a(T) — at, where 
t = (T — Tc)/Tc- By minimizing the Landau free energy 
Eq.0]with B = 0, we find the spontaneous magnetization 



M = \ — 



iit 



(5) 



Combining Eq. 0] and Eq. [5J the zero-field specific heat 
can be evaluated from the thermodynamic relation Cy = 
—T(d 2 F/dT 2 ), 



/n(Landau) 



a 2 T 
2bT 2 




forT<T c , 



for T>T C 



(6) 



Since this simple model only includes terms in the free 
energy that are dependent on the macroscopically or- 
dered moment, C m must vanish above Tc; the discon- 
tinuity in the specific heat at Tc is simply the limit- 
ing value as T -> T c , i.e. AC m = a 2 /2bT C - This 
value can easily be compared with experiment by not- 
ing that the gradient of the critical isotherm in the Ar- 
rott plot Fig. |21 gives 1/6 and that in the magnetiza- 
tion scaling plot Fig. 0] the interce pt o f the T < Tc 
data with the ordinate axis gives \J a/b directly. From 
the data we obtain b = 1.1 x 1CT 13 T (Am^ 1 )" 3 and 
\fa~Jb = 4.0 x 10 4 Am _1 . This gives a discontinuity 

AC'£ jandau ' ) = 150mJK -1 moL~ 1 , in excellent agreement 
with the experimental value of 155 ± 30 mJ K _1 mol -1 . 
The agreement between our measured AC and that cal- 
culated from the magnetization isotherms provides a con- 
sistency check on the Landau theory. 

The specific heat anomaly in weak itinerant ferro- 
magnets has been the subject of various theoretical 
studieo 53 i 54 i 55 ' 5 9i 5 7i 5 9. In the Stoner-Wohlfarth model 53 , 
the Landau coefficient a has the temperature dependence 



a(T) = -J± ( l - II 
{ ' 2 X0 V T 2 



(7) 



In this model the exchange field A is proportional to the 
magnetization; the T dependence of a reflects the reduc- 
tion in magnetization due to thermal spin-flipping exci- 
tations. The b coefficient can be obtained by minimizing 
F at T = 0, giving 



Mo 



where 



(dM\ 



(8) 



(9) 



TABLE II: Comparison of calculated and measured specific 
heat discontinuities. 





Landau 


Stoner 


SCR 


Exp 


AC(mJK~ 1 mol" 1 ) 


0.15 


0.33 


0.081 


0.155 



and Mo is the T = magnetization. As in the Lan- 
dau model, we can estimate the specific heat from F(T) 
[Ref.lH, 



Cl st ° 



Ml 



= 



T T 3 \ 

3^- forT < T c , 

T c If, J 

for T > T c .(10) 



Using M = 3.1 x 10 4 Am- 1 (= 0.17^ B f-u." 1 ) and 
Xo = 4.2 x 10~ 3 , as determined from M(H) measure- 
ments at T = 1.8 K [Ref. El, we find ACi stoncr) = 
330 mJK -1 mol -1 which is a factor 2 larger than the ex- 
perimental value. 

The Stoner-Wohlfarth approach can only include the 
effect of spin fluctuations through the renormalization 
of its phenomenological parameters. Experiments have 
shown that many of the properties of weakly ferromag- 
netic materials such as ZrZn2 and NI3AI cannot be ex- 
plained within this framework. Perhaps the most obvious 
property not explained by a mean field approach is the 
temperature dependence of the susceptibility above Tc, 
for which experiment shows a Curie- Weiss dependence 
X oc (T - Tc)" 1 as opposed to the x oc (T - T c y 2 
predicted by Stoner theory. In order to address the de- 
ficiencies of mean field theories, various self-consistent 
rcnormalized (SCR) spin fluctuation models have been 
proposed 3 * 5 ^. These follow the Landau- Ginzburg ap- 
proach and treat the local magnetization m(r) as a fluc- 
tuating stochastic variable. The SCR theory allows both 
magnetic corrections above Tc and the renormalizing ef- 
fect of spin fluctuations on the Landau a parameter to be 
taken into account. The effect of including spin fluctua- 
tions is to reduce the discontinuity in the specific heat 
AC m at Tc- This has been estimated by Mohn and 
Hilscher 55 based on the model of Murata and Doniach 5 , 



AC (SCR) 



MqMq 2 
4xo ?c ' 



Note that the Mohn and Hilscher value Ad 



(SCR) 



(11) 



=(1/4) 



ACm ton • The experimental value of AC m and the var- 
ious model predictions are summarized in Table ITT1 

In summary, the Stoner-Wohlfarth model overesti- 
mates the specific heat jump and the SCR-theory, as im- 
plemented in Refsi^ 5 -, underestimates the jump. It is 
pleasing to note however that the shape of the anomaly 
predicted by Stoner theory is similar to that observed ex- 
perimentally. Because of the nature of the theories one 
cannot attribute their failure to a single approximation. 
In the case of the Stoner theory the phenomenological pa- 
rameters (Tc, X0j b) used as inputs to the model are not 
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taken directly from band structure calculations, rather 
they are determined from the experimental magnetiza- 
tion M(B,T) and thus can be renormalized by fluctua- 
tions. The fact that the Stoner-Wohlfarth theory over- 
estimates the specific heat jump suggests that paramag- 
netic correlations give a significant contribution to the 
specific heat above Tc- The SCR-theory might be im- 
proved by using a more realistic model for the excitations 
(x"{q,u)) below Tc, where the longitudinal and trans- 
verse excitations are treated separately. Unfortunately, 
a full implementation of the SCR theory including, for 
example, the effects of the changes in the magnetic ex- 
citation spectrum on entering the ferromagnetic state is 
difficult. 



E. Conclusion 

In conclusion we have studied various thermodynamic 
and transport properties of the weak itinerant ferromag- 
net ZrZn2. Magnetization measurements show that the 



easy crystallographic axis is [111] and scaling plots reveal 
mean-field exponents for the temperature range stud- 
ied. Specific heat measurements reveal an anomaly whose 
shape is reminiscent of mean-field behavior. However, 
the measured discontinuity is significantly smaller than 
that predicted by Stoner-Wohlfarth theory. The resistiv- 
ity shows an anomaly at Tc and a non-Fermi liquid T 5//3 
behavior at low temperatures. Our results demonstrate 
the importance of collective spin fluctuations in ZrZn2, a 
material that was once considered to be a candidate for 
a Stoner ferromagnet. 
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